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The well-known D- dimensional Feynman integrals were shown, by Halliday and Ricotta, to be 
capable of undergoing analytic continuation into the domain of negative values for the dimension 
of space-time. Furthermore, this could be identified with Grassmannian integration in positive di- 
mensions. From this possibility follows the concept of negative dimensional integration for loop 
integrals in field theories. Using this technique, we evaluate three two-loop three-point scalar inte- 
grals, with five and six massless propagators, with specific external kinematic configurations (two 
legs on-shell), and four three- loop two-point scalar integrals. These results are given for arbitrary 

ON ■ exponents of propagators and dimension, in Euclidean space, and the particular cases compared to 

ON ' results published in the literature. 
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I. INTRODUCTION. 



The study of quantum field theories in the perturbative regime can sometimes become a very challenging and 
arduous issue to deal with, especially when one needs to go on calculating multi-loop Feynman integrals. There are 
several techniques that have been developed and refined in the course of time, but usually the calculations are rather 
involved and analytic results are difficult to obtain. Depending on the approach one uses to tackle the problem, the 
obstacles to surmount can be more or less demanding. For example, if one uses the standard Feynman parametrization, 
the integral over momentum space is straightforward — in the sense that nowadays one can easily fint them out even 
f^S " tabulated in textbooks on quantum field theories — but the resulting parametric integrals that are left over become 
ON [ quite complicated and cumbersome to tackle, especially with increasing number of internal momenta in the loops. Of 
course, if one is clever enough, he/she may be able to overcome some of these hurdles and even perform really hard 
| calculations 

There are, of course, other techniques which have been considered in the literature, e.g., integration by parts and 
Gegenbauer polynomial method in configuration space ||; the Mellin-Barnes representation for massive propagators 
D [ jl0| to cite just a few of them M], each one of those with its own strengths and weaknesses. 

• On the other hand, there is this novel technique known as negative dimensional integration method (ND1M), 
first devised by Halliday and Ricotta ||. It has as its starting point the principle of analytic continuation, it has the 
remarkable property of being equivalent to Grassmannian integration in (positive) D-dimensional space , it requires 
only the integration of Gaussian-like integrals and solving of linear algebraic equations and is suitable to deal with 
massless (|| or massive [|| diagrams in on- and off-shell regimes (?]] and even to carry out light-cone gauge integrals, 
with added troublesome gauge-dependent poles (|]. 

Our aim in this work is to demonstrate the simplicity of NDIM calculating some scalar three-point integrals, with 
five and six massless propagators, at two-loop level and some two-point three-loop integrals also in the massless case. 
The outline for our work is as follows: in section II we present a detailed calculation of a two-loop integral in the 
NDIM approach; in section III we write down the results for the remaining two- and three-loop scalar integrals, while 
in section IV we give our concluding remarks about this work. 

II. TWO-LOOP THREE-POINT VERTEX WITH SIX MASSLESS PROPAGATORS. 

This computation is performed following just a few simple steps Fist of all, let us consider the Gaussian-like 
integral, 

F = J j d D r d D q exp [~aq 2 - [3r 2 - j{q - rf - S(p + q) 2 - uo{q - k) 2 -6(q-r- k) 2 } 
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where the arguments in the exponential function of the integrand correspond to propagators in the diagram of Fig. 1 , 
with k 2 = t 2 = satisfying the on-shell condition. For the sake of compactness, we define 



A = a/3 + aj + a9 + (3-f + (35 + jS + 95 + (3w + •yuu + 9ui + (39. 

The parameters (a, (3, 7, 5, lu, 9) are quite general and arbitrary except that their real parts must be positive to make 
sure we have well-defined objects over the whole integration region. 

Expanding the exponential in the second line of (Q) in Taylor series and using the multinomial expansion for the 
argument of the exponential, we obtain, 

D ~ (-p 2 ) s *(-S:e - D/2)\a a /3 b 1 c 6 d Lj°0f 
F ~* , Xl l--- X4 l yi l---y 9 l Zl l Z2 l t-v-WW ( 2 ) 

where we have a fifteen-fold sum, with En = x\ + x 2 + x 3 + x 4 and we define, 

a = x\ + x 2 + xa + yi + y 2 + y 3 , 

b = xi + x 3 + yi + 2/4 + y 5 + y$ + z 2 , 

c = x 2 + x 3 + y 2 + J/4 + y 6 + y 9 , 

d = x\ + x 2 + x 3 + Xi + y 5 + y 6 + y 7 , 

e = ys + + z\ , 

/ = x 4 + y 3 + 2/7 + zi + z 2 . (3) 

On the other hand, taking the Taylor expansion of the exponential that is under the integration sign in (|l|), we 
have, 

°° ( _l\i+j+l+m+n+s 

F = E 1 ■ , , , oW^g'SW, (4) 
* — ' i\]U\m\n\s\ 

where S\j DIM is the relevant integral in negative D, defined by 

SUlM = Jd D q d D r(q 2 y(r 2 y( q ~r) 2l ( q +p) 2m (q-k) 2n ( q -r-k) 2s . (5) 

Now comparing both expressions for the original integral F we are led to the conclusion that, in order to have 
equality between these two expressions, the factor S^ DIM must be given by the multiple series, 

00 1 

S NDIM = p i \] — ; ; — ; ; — ; — T § a,i Sb,j <>c,l 5 d ,m <5 e ,« 5 f a , (6) 

^ xi\ ■ ■ ■ Xi'yi' ■ ■ ■ y 9 \zi\z 2 \ 

x,y,z—0 

with 

Pi = (-7r) D (p 2 ) CT r(l + »)r(l + j)T(l + l)T(l + m)r(l + n)T(l + s)T(l -a- -D), 

and for convenience we use the definition a = i + j + l + m + n + s + D. The system we must solve has, therefore, 
fifteen "unknowns" (the indices x, y, z) and only seven equations — six of which comes from comparing the exponents 
of (a, /?,■••) and one from the multinomial expansion, namely, Sy + Ez = — Ea; — D/2. So, such a system can only be 
solved in terms of some of these "unknowns" , and we conclude, from the combinatorics, that there are C>p = 6435 
systems (7 x 7) of linear algebraic equations, i.e., there are 6435 different ways of choosing these "unknowns". Of 
this grand total possible solutions, 4113 yield trivial solutions and for this reason present no interest at all. Yet, the 
remaining 2322 solutions are quite a big number! Fortunately, not all the remnant non-vanishing solutions need to 
be considered in order to solve the Feynman integral. In fact, only linearly independent series will be needed. In 
our massless case and where there is only one non- vanishing external momentum, the argument of such series is unit. 
Since NDIM provides, in general, hypergeometric series (see, for example ref. |8|] where we have double and triple 
series of such), we obtain in this case p F q (- ■ ■ |1), and from the theory of generalized hypergeometric functions p3| , ^4[ , 
a linear combination of linearly independent series of this type has up to p series. 

Consider then, the solution where the indices j/2, J/3, 2/4, J/5, J/7, ys , V9,Z2 are left undetermined. Solving this system 
we get, 
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xi = -I + m + n - s + y 2 + 2/3 + 2/4 - 2/5 - 2/8 + z 2 , 

x 2 =i + l + n + D/2 - y 2 - yg + y 5 + 2/7 + z 2 , 

x 3 =j + l + s + D/2-y 4 -z 2 , 

Xi = -n + s - y s - y 7 + y 8 + y 9 - z 2 , 

2/i = -m -n- D/2 - y 2 - y 3 ~ 2/4 - z 2 , 

2/6 = -o + m - 2/5 - 2/7 , 

z\ = n - j/8 - 2/g , (7) 

which substituted in (|^) yields an eight-fold summation. Seven of these can be rewritten in terms of gamma functions. 
The procedure is as follows: Firstly, the factorials must be converted into gamma functions and these into Pochhammer 
symbols, defined by, 

r \u\ - f \ r(n + fc) 

(n\k) = (n) k = , : ■ (8 
r(n) 

Secondly, use the properties, 

(n\j + k) = (n + j\k)(n\j), („| _ fe) = J=^_, (9) 

(1 — n\k) 

which follow from the definition above. Then, it is possible to rewrite some series in terms of Gaussian hypergeometric 
functions, which can, in general, be summed when its argument is unit Jl^ ]. 

vi h i-i a r(c)r(c-a-b) 

^i(M, c |i) = r(c _ o)r(c _ 6) - (10) 

Exception for this occurs when the parameters are such that c — a — b = 0. 

Therefore, if we rearrange these series conveniently, we can sum them, and applying this procedure for the series 
defined by indices 2/4, 2/3, 2/2, 2/8 > 2/9 , 2/7; 2/5; they all can be summed out. The last series, defined by index z 2 , however 
cannot be summed in this fashion since it is a 3-F2 function, which is summable only when it is of Saalschutzian type 
|jl3f , and this is not our present case. 

Many of the gamma functions that appear in this process simplify among themselves and in the end, what remains 
are fourteen gamma functions distributed in a fraction, i.e., seven in the numerator and seven in the denominator 
multiplying the z 2 series. This is the result in the negative dimension region and positive exponents (i,j,l,m,n,s), 
see ([|). Now we must be able to bring this result back, that is, analytically continue it, to our real physical world, 
that of positive D. This is carried out very easily: group the gamma functions in convenient Pochhammer symbols 
and use the second property 

Considering now the solution where Xi, X2 , £3 ,2/1,2/5, 2/6:2/9)^2 are left undetermined, we get another 
3-F2 (a, b, c; e, f\z) function. Let us recall that such functions have three singular points, namely (0,1, 00). Here, 
we are interested in solutions in the vicinity of z = 1. N^rlund Jl5[ proved that when the difference between e and / 
is an integer, some of the series concide or are without meaning. This is precisely our case. This means that our final 
result for the Feynman integral shall exhibit only two linearly independent generalized hypergeometric functions of 
the type 3 F 2 , instead of the usual three. We get, as a result, then 

Sndim = * D (/) CT {A 3^2 (01, 61, ci; ei, /1 1 1) + B 3^2 K , b[ , c[; e[, f{ |1)} (11) 



where 



and 



A = (-m\a)(a + D/2\ -2a- D/2)(-j\j + I + s + D/2)(-l\j + I + s + D/2) 

x {-i-j-l-s-D/2\i)(j + l + s + D\ - j - I - s - D/2 + m + n), (12) 



B = {-i\o){-m\a)(o + D/2\ - 2a- D/2) 

x - I - s - D/2)(j +l + s + D\-j- D/2)(-l - s\j + 1 + S + D/2), (13) 

where the set of parameters are defined as (ai, b\, ci; ei, fi) — (— s, m + n + D/2, —j — I — s — D/2; —i—j — l — s — 
D/2, l-j-s-D/2), and (a[, b[, c[; e' v f[) = (-a, -I, -j-l-s-D/2; -l-s, l + i-a). Note that {e x - fx) is 
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an integer. Observe also that, in the special case where all the exponents are equal to minus one, the hypergeometric 
function 3.F2 reduces to 2 Fi, which can be summed again, using (|Io|). So, for this special case of negative unity 
exponents, the final result is expressed in terms of gamma functions only. 

Let us point out here that our final result agrees with that given by Davydychev and Osland |llf| , i.e., meaning 
that equation (27) in r ef. § disa grees with ours. The difference is in the 9£(4) factor in which should correctly be 
read 9C(2). This gives the correct result for i = j = I = m = n = s = —1 and D = 4 — 2e. 

III. RESULTS FOR SOME TWO AND THREE-LOOP MASSLESS SCALAR INTEGRALS. 

The scalar integrals for the remaining diagrams (Fig. 2-7) are, in negative D, given by (number superscripts 2-7 
label the corresponding diagrams): 

°NDIM 



d D qd D r (r 2 Y(q 2 y{p + r) 2l (q + k) 2m (p + r-q) 2n , (14) 

4w« = J d°q d D r (r 2 Y(q 2 y(k q r f (k r) 2m (r+p) 2n , (15) 

SUlM = J d D q d D r d D k (r 2 Y(k 2 y(q 2 ) l (p - q - r) 2m (r - k) 2 \ (16) 

S% DIM = Jd D qd D rd D k(q 2 y(k 2 y(r 2 ) l (p-rj 2m (r-q) 2n (p-r-kj 2s , (17) 

S%DIM = I d°q d D r d D k (r 2 Y(k 2 y(q 2 ) l (r + k) 2m (p -k- r) 2n (r + k- q) 2s , (18) 



and 



SjfDiM = / d°k d D q d D r (q 2 y(k 2 y(p - k - r) 2l (k - q) 2m (p - k) 2n (r 2 ) s . (19) 



All these can be evaluated in a similar way following the steps of the previous section. For this reason we quote 
here only the final results, after continuation to positive dimension and negative exponents of propagators: 

SIdim = K D (k 2 ) p {C 3 F 2 (a 2l b 2 , c 2 ; e 2) / 2 |1) + D 3 F 2 (a' 2l b' 2 , c' 2 ; e' 2 , f 2 \l)}, (20) 



where 



and 



C = {-i\p){p + D/2\ -2p- D/2)(i + D/2 - p\p){p D/2) 

x -m-n- D/2)(-n\j + m + n + D/2) (21) 



D = (-i\p)(-l\p)(p + D/2\ -2p-D/2)(j + m + n + D\ -j-m-D/2) 

x (-j -m\-n- D/2){-n\j + m + n + D/2) , (22) 

with the set of parameters defined by (a 2 , b 2l c 2 ; e 2 , f 2 ) = (— m, i+D/2, —j—m—n—D/2; i+D/2 — p. \ — m—n—D/2) 
and (a' 2 , b' 2 , c' 2 ; e' 2 , f 2 ) = (—j, —p, —j — m — n— D/2; 1 + 1 — p, — j — m) and p = a — s. Note that here also we have a 
linear combination of two 3-F2 hypergeometric functions, since (e 2 — f 2 ) is an integer and N^rlund's theorem applies. 

Sndim = Tr D (k 2 y(-l\p)(p + D/2\ -2p- D/2)(i + 1 + D/2- p\p)(-j\j + n + D/2) 

x(-n\j + n + D/2)(j+n + D\ - 2j - 2n - 3D/2), (23) 



4 



&NDIM = K 3D / 2 (p 2 r'(-j\j +n + D/2)(-n\j + n + D/2)(-l\l + m + D/2) 
x (-m\l +m + D/2){p' + D/2\ - 2p' - D/2)(j + n + D\i) 

x(-j-n-D/2\-i), (24) 



where p' = p + D/2. 



Sndim = k 3D/2 (pY HI* + n + D/2)(-n\i + n + D/2)(-j\j + a + D/2) 
x(-s\j + s + D/2)(-i -1-71- D/2\l)(-j -m-s- D/2\m) 

x(cr' + D/2\ - 2d' - D/2){i + n + D\l){j + s + D\m), (25) 

where <j' = a + D/2, 

S%dim = ^ 3D/ Vr'HN + J + D/2)(-j\i + j + D/2)(-l\l + s + D/2)(-s\l + s + D/2) 
x(cr' + D/2\ - 2cr' - D/2)(-i-j -I -m-s- D\l + m + s + D/2) 
x(i + j + D\l + m + s + D/2){-n\ - I- s- D/2 + n) 

x(l + s + D\-l-s-D/2 + n). (26) 

Note that this result contains a special case (for m = 0) calculated numerically in Jig , Our result, on the other 
hand, was analytically obtained and is more general, since the exponents of propagators are arbitrary (cf. results in 
the second reference in ||). 

Finally, 

S 7 ndim = ^ D,2 {p 2 Y' + m + D/2)(-m\i + m + D/2){-l\l + s + D/2) (27) 
x(cr' + D/2\ - 2a' - D/2)(-i - j - m - D/2\j)(l + s + D\ - 21 - 2s - 3D/2) 
x(-l - n - s - D/2\a')(i + m + D\ - i + 1 - m + n + s)(-s\l + s + D/2). 



IV. CONCLUSION. 



Although NDIM is not a regularization method, it shares all the concepts of dimensional regularization since it is 
based on the principle of analytic continuation in the dimension of space-time. In this work we showed that with an 
adequate technique such as NDIM it is possible to calculate some two- and three-loop Feynman integrals with relative 
easiness and without much elaborate machinery, just the performance of some Gaussian-type integrals plus resolution 
of systems of linear algebraic equations. Of course, in view of the great number of systems of linear algebraic equations 
that are involved, we need the help of computer facilities, but the task is simple, which a PC can handle properly. 

The results we obtained can be used to study on-shell form factors in massless QCD at the two-loop level. Moreover, 
to illustrate how NDIM works beyond two-loops, we also performed some calculations of scalar three-loop, two-point 
integrals with arbitrary exponents for propagators and dimension, in Euclidean space. 



ACKNOWLEDGMENTS 



A.G.M.S. gratefully acknowledges FAPESP (Fundagao de Amparo a Pesquisa do Estado de Sao Paulo, Brasil) for 
financial support. 



[1] G.Kramer, B.Lampe, J.Math.Phys. 28 (1987) 945. 

[2] R.J.Gonsalves, Phys.Rev. D 28 (1983) 1542. W.L. van Neerven, Nucl.Phys. B 268 (1986) 453. 
[3] I.G.Halliday, R.M.Ricotta, Phys.Lett. B 193 (1987) 241. 



5 



G.V.Dunne, I.G.Halliday, Phys.Lett. B 193 (1987) 247. 

A.T.Suzuki, A. G.M.Schmidt, JHEP 09 (1997) 002. A.T.Suzuki, A. G.M.Schmidt, Eur.Phys.J. C5 (1998) 175. A.T.Suzuki, 
A. G.M.Schmidt, Phys.Rev. D58 (1998) 047701. 
A.T.Suzuki, A. G.M.Schmidt, J.Phys.A31 (19 98) 8023. 



A.T.Suzuki, A.G. M.Schmidt, hep-th/9712104 



A.T.Suzuki, A. G.M.Schmidt, R.Bentm, Nucl.Phys.B537 (1999) 549. 

K.G.Chctyrkin, A.L.Kataev, F.V.Tkachov, Nucl. Phys. B174 (1980) 345. K.G.Chetyrkin, F.V.Tkachov, Nucl.Phys. B192 
(1981) 159. 

E.E.Boos, A.I.Davydychev, Theor. Math. Phys. 89 (1991) 1052. A.I.Davydychev, J. Math. Phys. 32 (1991) 1052; J. Math. 
Phys. 33 (1991) 358. 

A.I.Davydychev, POsland, Phys.Rev.D59 (1999) 014006. 

A.Erdelyi, W.Magnus, F.Oberhettinger and F.G.Tricomi, Higher Transcendental Functions (McGraw-Hill, 1953). 
I.S.Gradshteyn, I.M.Rhyzik, Table of Integrals, Series and Products (Academic Press, 1994). 
W.N.Bailey, Generalized Hypergeometric Series (Cambridge Univ. Press, 1935). 

Y.L.Luke, Special Functions and Their Approximations (Vol.1, Academic Press, 1969). L.J. Slater, Generalized Hypergeo- 
metric Functions, (Cambridge Univ. Press, 1966) 
N.E.N^rlund, Acta Math. 94 (19 55) 289. 



R.Easther, G. Guralnik, S.Hahn, hep-ph/9903255 



p 





1 



r + k 




p — r — k 



p — k — 



Figure 7 



1 



